Electronic phase transitions in the half-filled ionic Hubbard model 
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A detailed study of electronic phase transitions in the ionic Hubbard model at half filling is 
presented. Within the dynamical mean field approximation a series of transitions from the band 
insulator via a metallic state to a Mott-Hubbard insulating phase is found at intermediate values of 
the one-body potential A with increasing the Coulomb interaction U. We obtain a critical region in 
which the metallic phase disappears and a novel coexistence phase between the band and the Mott 
insulating state sets in. Our results are consistent with those obtained at low dimensions, thus they 
provide a concrete description for the charge degrees of freedom of the ionic Hubbard model. 
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One important problem in the field of electronic struc- 
ture theory is the understanding of the Mott metal- 
insulator transition (MIT)jii 2 - Many details have already 
been clarified using the canonical model for this tran- 
sition — the one-band Hubbard model. At half filling 
(i.e., having in average one electron at each lattice site) 
this transition from the metallic to the Mott-Hubbard 
insulating (MI) phase occurs with increasing the on-site 
repulsion U. Great progress in understanding the MIT in 
Hubbard-like models has been achieved in the last decade 
with the development of the dynamical mean field theory 
(DMFT) . Although many aspects of the coherent (Fermi 
liquid) metallic phase and the incoherent MI phase are 
now quite well understood, other questions remain open, 
especially the relationship between the MI regime and 
band insulators (BI), where one sub-band is almost com- 
pletely filled by two electrons such that an excitation gap 
is formed. 

To study the interplay between band and Mott- 
Hubbard insulators, various extended versions of the one- 
band Hubbard model have been proposed. It is remark- 
able that different models show separate kinds of behav- 
ior. Some of them have a crossover from the BI to the MI 
regime,- whereas others have clear transition points with 
a metallic phase in between^ A widely studied model for 
the second class is the ionic Hubbard model (IHM)£&£ 
In the one-dimensional (ID) IHM, however, a ferroelec- 
tric (or bond-ordered) phase is realized which separates 
BI and MI, and the metallic phase shrinks to only one 
point i 10 ' 11 Already in ID, a finite metallic region can be 
recovered by introducing an intra-sublattice hopping t' 
into the IHM»i2 In 2Di or at larger dimensions 5 a cor- 
relation induced metallic phase was reported between BI 
and MI, but it is still under debate whether this metallic 
phase shrinks to a line 5 or if it ends up at certain par- 
ticular point by increasing the ionicity A. In this work 
we show that the metallic phase disappears at a certain 
value of A above which we found a coexistence region 
composed of distinct insulating phases. 



It might be not very surprising to find quite different 
behavior in various models at different dimensions. But 
even if we restrict our attention to the IHM, the phase 
diagram is highly disputed. The 2D case was studied 
by a cluster extension of DMFT— and by the determi- 
nant quantum Monte Carlo method^ showing different 
results. The cluster-DMFT study suggests a bond-order 
phase separating BI and MI regimes, while the stabil- 
ity of the metallic phase was established on a 2D square 
lattice at finite temperatures." 

To clarify the many open questions related to the elec- 
tronic phase diagram of the IHM, we study this prob- 
lem here in the limit of high dimensions. In this limit 
the self-energy is k-independent and, thus, the dynami- 
cal many-body problem can be mapped onto a single-site 
one. Using basically the same DMFT scheme as Garg et 
air- we confirm the existence of the metallic phase for 
moderate values of U/A, but we find also important new 
features which are different from those reported in this 
recent work. So, for intermediate values of A, we find 
a continuous Bl-to-metal phase transition and a discon- 
tinuous one between the metallic and the MI phase with 
increasing the Coulomb interaction U. 

We start our study with the IHM 

H= -t (ctc^+H.c.) + Uj2^nn 
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+ Aj^h,- Aj^h, (1) 
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on a bipartite lattice with two sublattices A and B hav- 
ing different on-site energies ea = +A and eb — — A. t is 
the (inter-sublattice) hopping term and U is the on-site 
Coulomb repulsion. cf a (respectively Ci a ) denotes the cre- 
ation (respectively annihilation) operator of an electron 
at the lattice site i with spin a and rij ff is the occupation 
number operator for electrons with spin a . We consider 
an average number of one particle per lattice site, i.e., 
the half filled case. 
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Application of the DMFT scheme to the IHM implies 
one to solve two (A, B) local problems given by the local 
one-particle Green's functions G a . At high dimensions, 
these local propagators are connected self-consistently by 
the Dyson equation = G^ 1 + S Q . 9 Without loss of 
generality in what follows we will perform our calcula- 
tions using a semi-circle density of states (DOS) — this 
type of one-particle DOS is realized on a Bethe lattice 
with large coordination number. 

On a bipartite A—B lattice with nearest- neighbor hop- 
ping t the single-site Green's functions reads 



G q (cj) = £ a 



°° ale- 



(2) 



with a — A (a — B) and £ Q = u) + i0 + — e a + 
/i — T, a (uj). After integrating Eq. ((2]) using the semi- 
circular DOS po(e) = V4t 2 — e 2 /(2irt 2 ), we obtain two 
coupled equations Go,(w) _1 = £ a — t 2 G a (w) for the IHM 
on a Bethe lattice. As in Ref. @, we use the iterated 
perturbation theory (IPT)i£ to calculate the local self- 
energies E Q (w). From a direct comparison between dif- 
ferent methods£ ' 13 i 14 ' 15 it is found for the one-band Hub- 
bard model that the quasi-particle peak (Kondo reso- 
nance), the position of the Hubbard bands, the Mott- 
Hubbard gap, as well as the principal characteristics 
of the metal-insulator transition in infinite dimensions 
are all well captured by the IPT solver. Since at each 
lattice site the band filling is always away from one- 
half in the IHM (see Fig. 0J, it is proper to use the 
canonical IPT ansatz: = + A^a {uj) / {\ - 

B a Ti a (w)) for the correlated problem. Within this 
ansatz E^f F = Un a is the first-order Hartree term 
and n a = —1/ttJ_ lmG a (u>)duj. The coefficients A a 
and B a are determined by requiring S q (cj) to be ex- 
act at high-frequencies and in the atomic limit. Finally, 



Y,a\uj) = S„ '(uj, [Qo a {^>)]) is the well known second- 
order diagram for the self-energy^ 

In our treatment, the host Green's functions Q$ a re- 

( 2^ 

quired to compute S„ (ui) self-consistently are written 
as Qoa = Goa + £q — The energy shift e a — [i intro- 
duced here accounts for the correct atomic limit of the 
IHM at U ^ A, t (see below) as well as the mirror-type 
particle- hole symmetry between sites A and B. Notice 
that a similar type of correction is employed in the con- 
text of single-site Hubbard model (ea = £b = 0), where 
the energy shift is — \i (= —U/2 at half filling)^ In addi- 
tion, we point out that our energy shift is different from 
the Hartree correction — £^f F employed in Ref. HI, and 
this might be the cause for significant differences between 
our and their results in the large U limit. 

Let us now present our results. In Fig. [TJ we show the 
spectral density for the charge phases of the IHM. In the 
upper panel, we present our result for the BI phase for 
U/t = 0.1. The full line corresponds to the two-site av- 
eraged (or the total) DOS, whereas the dotted line is the 
site resolved B-DOS. As expected the charge gap is small, 
of the order of A. We observe clear van-Hove singularities 
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FIG. 1: Total and site selective (B) density of states (DOS) 
for A = O.lt and various values U. The parameters are chosen 
identical to those of Garg et a/.,— leading to very close results 
besides for the lowest panel (larger gap and V-shape leading 
edge in the DOS). Notice the small difference between the 
lower and the upper Hubbard band in the site-selective DOS, 
this results from the softening of the CDW at large U. 



at the lower band edge at each a-site. With increasing 
U the band gap is suppressed (by dynamical transfer of 
spectral weight characteristic of correlated electron sys- 
tems) and a metallic phase sets in. As can be seen in 
the dotted line of Fig. [l]-(b), the B-DOS shows a pseudo 
gap feature at the Fermi level. This implies that the one- 
particle DOS in the metallic phase is not pinned to its 
unperturbed (U = A = 0) valued. This is understood, 
because the charge density wave (CDW) regime prevents 
the Fermi liquid fixed point at any value of the Coulomb 
interaction. As expected in such cases, the imaginary 
part of S a (w) (not shown) is finite at the Fermi level. 
The dynamical transfer of spectral weight found in the 
metallic phase can be traced to a Kondo like process that 
leads to strong electron mass enhancement. A similar ef- 
fect is known to take place in correlated electron systems 
like the periodic Anderson model^ or in some regimes of 
the disordered Hubbard models By further increasing 
U, the Mott-Hubbard insulator with a gap of the order 
of U — W (W is the bandwidth) sets in as the third phase 
of the IHM. 

Having chosen exactly the same parameters as in Ref. [H 
(and basically the same method of solution), we find very 
similar results for the upper two panels, but substantial 
differences for the results in the large U regime are ob- 
served. In particular, the van-Hove structures found in 
Ref.! near the Fermi energy broadens and shifts to high 
energies for U — 7.5i. There is an appreciable modifi- 
cation of the whole lineshape. Within our self-consistent 
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FIG. 2: Charge gap of the IMH as function of [/ for various 
values of A. The coexistence region between metal and Mott 
insulator (MI) is reminiscent to DMFT for the pure Hubbard 
model.— At large values of A the metallic phase shrinks and, 
we obtain a coexistence region for the band (BI) and Mott 
insulators. 



scheme the Hubbard bands centered at w w ±4t are in 
very good agreement with those found for the one-band 
Hubbard model. Clearly, our results for U = 7.5t are 
more plausible than those obtained in Ref. d. Moreover, 
it is understood that in the large U limit the staggered 
charge density is small (see results of Fig. 2]), forcing one 
to recover to a good extent the physics of the regular 
Hubbard model. This is exactly what is shown in the 
lower panel of Fig. [TJ 

To characterize the correlated nature of the various 
phase transitions and the stability of the metallic phase, 
we show in Fig.[5]the obtained charge gap as a function of 
U. As can be seen in this figure, the width of the metallic 
phase shrinks until it disappears at a critical value. For 
A < 1.9i, we find a continuous transition from the BI to 
a unique metallic phase with increasing U. Differently 
from the Bl-to-metal transition, the transition from the 
metallic phase to the Mott insulator phase is of the first 
order. The latter transition in the IHM is completely 
identical to the corresponding transition in the ordinary 
Hubbard model. Across the metal-to-MI transition line 
we observe a coexistence region consisting of two distinct 
(numerical) solutions, which is accompanied by a hystere- 
sis loop characteristic of a first order MITfS In agreement 




FIG. 3: Phase diagram of the IHM obtained by single-site 
DMFT. The dashed area corresponds to the coexistence re- 
gion between metal (M) and Mott insulator (MI) or between 
band insulator (BI) and Mott insulator, respectively. The 
line separating the BI and M phase is characterized by the 
vanishing of the charge gap. 



with previous theoretical studies^ the width of the co- 
existence region decreases with increasing ionicity (A). 
For larger values of A, the metallic phase disappears and 
a coexistence region between BI and MI is found. 

The obtained phase diagram for the IHM on a Bethe 
lattice is shown in Fig. [3] It differs considerably from 
those found and/or recently proposed for the IHM on a 
usual bipartite lattice in two or at high dimensions £&2& 
Differently from previous single-site DMFT results^ the 
metallic phase is now considerably widened up to 1.9t. 
The first-order metal-to-MI and BI-to-MI phase transi- 
tions resemble the coexistence region found by Kancharla 
et al. in 2Di£ In addition, our window for the metallic 
region at moderate values of A is similar to the deter- 
minant quantum Monte Carlo study also in 2D, 6,8 indi- 
cating common similarities between the physics at low 
and high dimensions for the IHM. To further illustrate 
the underlying features in the coexistence region and the 
first-order phase transitions, we show in Fig. [4] the stag- 
gered charge density ub — tia as a function of U . As one 
can see, there is no peculiarity at the transition point 
between the BI and the metallic phase. However, the 
staggered charge density jumps at the transition towards 
the MI phase. 

The possibility of having different types of first-order 
phase transitions in the IHM (see Fig. [3J is one relevant 
aspect of our work. Furthermore, our results show that 
the coexistence region is not necessarily connected with 
the bond-order phase recently found with cluster-DMFT 
for the 2D-IHM. 7 We find that a first-order transition 
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FIG. 4: Staggered charge density nj — tia as a function of 
(7 for different values of the ionic potential A. Our results 
for A = 2t are in good agreement with those obtained within 
cluster-DMFT in 2D.- 7 - 



between metal and MI (or between BI and MI) exists 
not only in 2D but also at high dimensions. Thus, the 
first-order MIT of the IHM is of the same nature as that 
found within the single-site DMFT for the one-band Hub- 
bard model. 9 Finally, we recall that the 2D 7,8 works reach 
opposite conclusions regarding the metallic phase of the 
IHM, indicating a great sensitivity on the procedure used. 
The metallic phase derived within determinant quantum 
Monte Carle- suggests a quasi-particle peak at ui = 0, 8 
while a pseudo-gap like feature is seen in the Lanczos so- 
lution of the cluster bathi 7 - Clearly, near the Fermi level 
our metallic phase is more consistent with the CDMFT 
calculations. However, the shape of our phase diagram is 



similar to the one derived in Ref. @ which indicates that 
the correct nature of the metallic phase in 2D remains to 
be investigated. 

From the pure Hubbard model it is known that IPT 
method is well suited to show the existence of the coex- 
istence region and provides a correct qualitative picture 
of the metal insulator transition in infinite dimensions. 
However, there are some details of this transition, like e.g. 
the precise position or the stability of different phases 9 -, 
which have to be clarified by more sophisticated methods 
and which are left for further studies. 

In summary, we derived the many-body phase diagram 
of the IHM in the limit of infinite dimensions. For inter- 
mediate values of the ionic potential A, we find a metallic 
phase between the BI and the Mott-Hubbard insulating 
phase which disappears, however, for larger A. Our re- 
sults clearly differ from a previous work which employs 
the same single-site approach for the IHM£ We believe 
to have found the correct spectral density in the limit of 
U 3> A. Differently from Ref. [5j, our result in this regime 
shows clear Mott-Hubbard features, in agreement with 
the general expectation*?- We find first-order phase tran- 
sitions between the metallic and the MI phases as well as 
between the band and the Mott insulator. The latter one 
were not reported in previous studies on metal-insulator 
transitions for the IHM. We show that a discontinuous 
transition exists already at high dimensions and is not 
necessarily connected with the bond-order phase as pro- 
posed in Ref. 0- 
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